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ment of Mathematics and Statistics, Bundelkhand University, 

Jhansi. 

This THESIS consists of three Chapters numbered I, II 
and III; every Chapter is divided into several sections 
(progressively numbered like 1.1, 1.2, 2.1, 3.1 ' • 

The formula and equations are numbered progressively within 

every section. 

Chapter I has been divided into several parts, first 
is History of Fluid Variables and others are the Fundamental 
Equations described in detail. The existence of surface of 
discontinuity and other properties of shock wave is detailed 
out. The generalised Rankine-Hugpniot Jump relation is derived 
in the preceding section. Similarity principle is defined sue 
as the partial differential equation reduces to ordinary diff- 
erential equation with independent variable. The structure of 
a radiative shock wave in a viscous compressible thin layer 
adjacent to the surface of body has been studied and expressions 
for the variation of the physical quantities across the shock 

have been derived. 

The present work is with well-known result notation 
given at the beginning of the Thesis. References are given 
at the end of the Thesis and are denoted in Thesis of numbers 


in square brackets® 

The two papers are solved in Chapter II and Chapter 
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(U.P. ) , for his valuable suggestions. 

My special thanks are to Dr. V.D.SHARMA. Department 
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CHAPTER-I 


INTRODUCTION 


Fluid Dynamics is that branch of Science which -is 
concerned with the study of motion of fluid or that of 
bodies in contact with liquids and gases. Which are classi- 
fied as fluid. We regard liquid as a incompressible fluids 
for all practical purpose and gases as compressible fluid. 
In general, fluids have five physical variable? Density, 


Velocity, temprature, pressure and velocity 


M f 2 ] 


Since the phenomenon considered in fluid dynamics 
are macroscopic, a fluid is regarded as a continuous medium# 
this means that any small volume element in the fluid is 
always supposed so large that it still contain a many number 
of moleculars so when we take of infinitely small element 
of volume. We always mean those which are physically small 
it is found to consist of molecules in random motion and 
separated from one another by distance which are at least 
comprable with molecular size. In the case of gases, the 
separation distance are great? in the case of liquid, they 
are less great and in the case of solids even less so. 

The mathematical description of state of a moving 

fluid is affected by means of function which gives the 

distribution of fluid velocity V=v(x,y, z,t) end. of any two - 

thermodynamics quantities pertaining to the fluid. For 

example the pressure P (x, y, z, t) and the density -f (x,y, z,t) 

All the thermodynamics quantities are determined by the 

* • , 

values ©f any two of them together with equation ®f state. 



Hence if we are gives five quantities namely the three 
components of velocity v, the pressure p, and the density f 
. The state of moving fluid is completely determined 

M h) 


There are two method fox finding motion 


mathematically, these are » LAGRANGIAN • and “ EULERM 


method for refer to « INDIVIDUAL TIME RATE OF CHANGE 


and • LOCAL TIME RATE OF CHANGE “ respectively. 


They are as 


(1) “INDIVIDUAL time rate of change" 


In this method we study the history of each particle 
i»e any fluid particle is studied and is pursued on its 
onward course. Observing the change in velocity, pressure, 
and density at each point and at each instant* 

The fundamental equations ©£ motion in Lagrangian 
form are non-linear and hence it leads to many difficulties 
■while solving a problem* In fact, it is employed with am 
advantage only in some one-dimensienal (involving one space 

coaftfitt&A-e) Pfohlr .. {'] [ c ] 


(2) “LOCAL TIME OF RATE OF CHANGE 81 


in this method we study a general point in space 
occupied by fluid and study the whatever changes taXe place 
in velocity, pressure, and density as the fluid posses 
through this point rather than the variation of velocities. 


and acceleration etc© 
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pressure 


As they follows their owe path® in it tn«s 
luid particles are not identified instead a point in chosen 
.nd changes in velocity etc. are studied as the fluid passes 
h rough the chosen fixed point® £ 5 J 


We can drive the 'LOCAL 8 and 8 INDIVIDUAL * time rates 
following or relationship between Lagrangian and Eulerian 


nethod. 


Let the fluid motion be associated by a scaler function 
-=fHx,y,z,t) or «=^(r, t) keeping the point P (x,y, 2 ) fixed. 

Che change in during an interval of time ' Sir* is 






° r ^ 

;ence the' local time rate of changes* is given by 

1, ~ 'Jzl 


Now keeping the particle fixed the change in <=^ is 

( ft it , ^ t&lr) — 

When is the change in position of fixed particle 

uring the short time * therei-ore# 






fib- 




Given the m individual time rate of change 


Let q (u, v,w) be the velocity of fluid particle# 
ach that 

~zr* a 
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~ r* ft 0 a . 

U 1 -f- ■'f 


and 
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etc. 
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This is relation between the two time rates. 
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FUND AM E NT AL EQUATIONS GOVERNING THE FLOW OF GAS 


We consider the fluid fl v .a wish to fxnd the 
velocity distribution as well as the states of the fluid 
over all spaces for all time. If the fluid is a single 
gas ordinarily a knowledge of the three velocity components 
(u, v,w) the temperature T, the pressure P # and the density 
of the fluid which are function of spatial coordinate and 
time is desired. Hence we study three dimensional flow 
and consider all three components of velocity together with 
the pressure, density and temprature of a fluid as a function 
of three spatial coordinates x, y and z and the time t. 

We shall consider only inviscid and non heat conducting 
fluid because the effects of viscosity and heat conduction 
are usually negligible except near a solid boundary ©r inside 
a shock® J 

(a) EQUATION OF STATE 


Which connected the temperature T, dot pressure P 
and the density of the fluid, for a perf act gas, this 
equation may be written as 


¥^ 1*7 


£'■ 'J 


R being a gas constant® 

(b) EQUATION OF CONTINUITY 

Which express the conservation of mass of the fluid. 
If q represent the velocity of the fluid at any time t* the 


£>2 7 


equation of countinuity may be written as 

JL. h- oL' v ( -P ? ) *d 

d>t" 

(c) EQUATION OF MOTION 

Which express the relation of conservation of momentum 
in the fluid. Neglecting the body force and considering only 
the inertial forces and pressure forces, the equation of motion 


:n by 

C> - 3 ] 


the fluid. 

tM.J 

Where e and Q are the internal energy and the energy 
generated by external sources per unit time per unit mass of 
the fluid. 

(e) MAGNETOHYDRO DYNAMICS EQUATION 

Under magnetohydro dynamics we study the motion of 
electrically conducting fluids in the presence of electro- 
magnetic fields. The equation of motion of fluid element 
from fluid mechanics with the macroscopic phenomena logical 


may be written as ^ 2 1 

i>V 


Where 


'P 


e»— oj Jo 


A- 


is the usual mobile operator < 


£> 


J>kr 


o 

c)4~ 


i~ C J 


(d) EQUATION OF ENERGY 

Which express the conservation of energy in 
This can be expressed as ^2 "J 


$Jr 
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f J- ) 

-t 

1 7>Jr 

L < J 
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maxwell equation from electro dynamics and the thermal 
ami # equation of state from thermodynamics. The eleetrie 
current is induced in a conductor moving in a magnetic 
field. These currents being in a magnetic field suffer a 
mechanical force called Lorentz force. 

Magnetohydro Dynamics its peculier interest and 
difficulty to this interaction between the field and the 
motion. 

Neglecting maxwell* s displacement currents and taking 
magnetic permeability Jx as unity. The field equation are 

Cuff H - -f-£ '( dP“l dt-J 


Cud £ 




V £ 


cX'v H 


_ 1 

— £ 
dV 

c 

d*r 

M < 
Gr 

. ^ 

% 
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~ — 7 ^ 

E 

and 


f i-*! 


magnetic field, current density, electric intensity and 
velocity of light respectively, £ is dielectric constant 
and the charge density. In addition the constitutive 

equation are __ 

TAtr r , ^ e 
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“ ORIGIN OF DISCONTINUITY " 

The laws of conservation of mass, momentum and energy 
that from the basis for the equation of inviscid flow of a 
non-conducting to not necessarily assume continuity of the 
flow variables. These laws are originally formulated in the 
form of differential equations. Simply because it is assumed 
at the beginning that the flow is continious. These laws, 
however can also be applied to these flow regions where the 
variable undergo a discontinious change, from the mathe- 
matical point of view, a discontinuity can be regarded as the 
limiting case of very large but finite gradients in the flow 
variables across a layer whose thickness tends to zero. 

Since the dynamics of an inviscid and non-conducting 
gas there are no characteristic lengths, the possibility of 
the existence of arbitrary thin transition layers is not 

l 

excluded. In the limit of vanishing thickness, these layers 
reduce to discontinuity such discontinuity are called SHOCK 

WAVE. 

Actually the shock wave is not a simply surface of 
discontinuity but a very narrow region in which a large vari- 
ation of pressure and velocity occurs it is possible to estimate 
the thickness of a shock wave approximately. From this appro- 
ximate estimate we may show that the extent of the shock thick- 
ness is very small and therefore neglected and may be consi- 
dered as a surface of discontinuity P t [ 2 1 t 4 l l 6 ] 

Shock wave appear in many comp ras sib le flow problem 

- wlth neat addition is the “CONDENSATION SHOCK 

or ten a sn^cx wave witn neac j 


In the actual flow of air, there is normally a certain 
amount of moisture present which is a vapour state. Second 
one is heat addition in "DETOKATIOL WAVE" which arises from 
the rapid transformation of explosive material and the flow 
direction is not perpendicular to the shock wave front such 
in the case of an "OBLIQUE SHOCK". There is a simple 
graphical relation for the velocity components across an 
oblique shock known as the shock Polar. 


t 
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THE MATHEMATICAL TREATMENT OF SHOCK WAVE: 




1st 

S ■ 

1 §1 
ri-f 
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Ws c_n define a density ^ unit length and 

a flux ^I'X/ 1) per unit time® We can define a flow velocity 


V" C - *; *) by 


V 


‘ilf 


We can stipulate that the rate of change of the 
total amount of it in any section ~y X J ^- 2 . must be 
balanced by the net inflow across % ( and • That is 


o(_ 

If 


01, 


?(- 


f> 0,t)cU + ^ (.*■, t ) - c \,L'x>" M ) — 0 


-f (•», *■) 


has continious derivatives, we may take 


the limit as ^ and obtain the conservation equation, j j 


il 


-h 


iL- 


The simplest wave problems arise when it is reasonable;' •’ 
on either theoritical or empirical grounds, to postulate a 
functional relation between q and ^ „ 

1 - d(f) [j 

Consider first the mathematical question of whether i 
discontinuities are possible. Suppose there is .a discontinui-' | 
at x = S (t) and that 'Xj and are chbsen so that f 

| and ^ and their first deri- 1 
"X , X > £ ( ir) and in S (Jr) >*-'>■ ^2 


E.i; 

5 

■i 

■■ 


H | > S CO>*»- suppose 
vatives are continious in 





mmf0i 


They will show that weak solution of a conservation law 


It is the true conservation equation then it may be 
deduced as the shock condition by the same argument. That 
correct choice of weak solution is made on the basis of 
which quantities are really conserved across the shock. 
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TYPES OF OBLIQUE SHOCK WAVE 

(1) Plane Shock 

(2) Cylinderical Shock 

(3) Spherical Shock 

One of the interesting case of unstudy an isentropic 
flow in the decay of a strong shock we may consider such 
a problem as depending on time t and on a single spatial 
coordinates, he. we may consider either a plane (or normal) , 
cylindrical, or spherical shock wave £ 4 ] ^6 J 


( 1 ) 


PLANE SHOCK 


Consider a normal shock to oe a surface of discontinuity 
in velocity, pressure, density and temperature of the fluid. 

We shall find the relation between these quantities in front 
of and behind the shock, we cnoose the coordinate system such 
that the shock is stationery and the fluid moves it. Under 
these conditions the flow is steady, we assume that the fluid 
is an ideal gas, so that the perfect gas law holds both infront 
and behind the shock, and the specific heat are constant, we 
assume that the velocities of fluid particle are perpendicular 
to the shock front. 


( 2 ) 


CYLINDRICAL SHOCK 


If there is an instantaneous release of energy along a 
line. The shock front such developed are known as cylindrical 


shock wave 


• ' j-jar t" w " 1 i 1 
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k 3) SPHERICAL SH nr r 

There is an instantaneous release of energy along 

1 POlnt • The Shock Trent such developed are known as 
pherical shock wave© 



LAW OF MOTION OF GAS PARTICLE 

A shock wave may be produced by an instantaneous 
:rgy release* such as the explosion of an atomic bomb, 
assume that a finite amount of energy is suddenly 
eased in an n rtf? nttely concentrated fora in the spherical 
e, and as energy per unit length E in the cylindrical 
«o If the value of energy release E are the same 
both cases© 


' ilf l. 
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We discuss here the fundamental equation in detail 
which gives the basic inpastructure to the problem elaborated 
in the later part of the thesis,, 


(1) EQUATION FOR THE LAW OF CONSERVATION OF MASS 

When a region of a fluid contains neither sources 
nor sinks that is to say when there are no inlets or out- 
lets through which fluid can enter or leave the region. 

The amount of fluid within the region in conserved in accor- 
dance with the principle of conservation of matter £ 1.3 



v "\ 


Let ZXS be a closed surface drawn in the fluid and 
taken fixed in space and let -I — (x, y, z, t) be the fluid 

density (mass per unit volume) at any point (x,y, z) of 
the fluid in at any time t. Suppose n is the unit 

out ward drawn normas at any surface element ^5 of & 5 
where ^ £ <£ &■$> , then if q is fluid velocity at the 

element 8S * The normal component of q measured out ward 
from is n.q. Thus 

Rate of efflux of fluid mass per unit time 

^ ^■%SS 


mt 


across 


mm 


M— 



ox mass riow out of £x 


At time t, the mans o£ fluid with in the element 


Local rate of mass increase within A.\> 


for all volume 




r 




“ - - 1 ■< ■ ■ . • r 
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tion C $3 s. the general equation of continuity 
sources and sinks. Since 

** Hi)- f v \ -t i v f 

Substitute the value oj = v/ / ad 1 • , 

6 ° f ^ <• « equation CrS) then 


il 

af 

K 

ar 


-+ 


•f • v i -t- V v / 


— 1> 


+• j y v =- 0 


5P 


H- V • ^ & 


LmJ 

£M.] 


In second and third equation d/rft- denots 

differentiation following the fluid motion and the 
operational equivalence 




A -t i ' v 


has been used. 


. the EQUATION OF MQT Tnw 

When the fluid moves an electric and magnetic field 
the body force F per unit volume oonsiste of three parts. 
Gravitational, Electric and Magnetic, the gravitational ' 

body force per unit volume is f f' where g is the accel- 
erati on cbje to <**??•'?' % n - * 

W * ' * * ' * **“* V Ox Uil&c ^ Q£ th$ 

fluid would c °ntain a charge of amount ^ S^’so that the 
force on it due to an electric field by intensity E would 


nora,ai cross section of a fluid element whose 
length $ S is the direction of j. This element moves 
along with the local fluid velocity v in a magnetic 
field of intensity H» 



kinetic coefficient of viscosity 

for ordinary viscous fluid 

^ — p- - £ -t v ip 

If the fluid is under the effect of magnetic field 

J)\A 




h a _ 


j,g - wijb + o Vi v 


: i ii _ 
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When external force -f Sc mechanical force -f ^ ' has 
been neglected. £> — 


r~> 


•f> J)\* 


(Jui A p ) - -t3v*r> 


As 


cj< 


tluyJ M •=- u + 


c3 sf ’ 


M < 


■ 

I 




neglecting the displacement current 
Curl H ^ 3 

Substitute it 

^ - t-Ti + ' 3 ' g ' v ' 

by vector identity f 

= C ** * o*r.< h* ) •+ ( h 5 ^ 7 1 ^ 

^ O^riXtP)^ 


Substitute it 


J . __q » . j^L 

^ ^ 'T r ‘^ 1 -71 + w* 


ilff’vW -t^'% 

,J ' g<P 


^ ^ vj (h %t) t. 4 

PA- 1 0 ' I 


PA— 

Where |p I- 


( b + Al ' lL 


pressure. 


3A 


is magnetic 


-P 


Pje 


c jr r<=< M *4, A £ iT° s? ) 

0 _i_ i7 ^7*2- 


^ is coefficient of viscosity 

This gives the equation of motion f®r viscous 
conducting fluid. 

(3) EQUATION OF ENERGY 


-t <7 






For one dimentional steady flow of an inviscid 
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compressible fluid in a channel, the momentum law gives 

<d f 


f 


% dl 


■'t> 


here q is the velocity, p is the pressure and "f is 
the density of the fluid. 

Since we consider only inviscid and non heat conducting 
fluid, the conservation of energy gives the energy equation as 

cl 8 Tills ^ 4e -f M t y o 
d8 j= dh -f' ^ P ' ~ 


which, in vector notation, becomes 




vH - j ^ \° 


it is sometimes convenient to introduce a stagnation 

enthalpy Ho defined as 


k 


H + 1 


L’"*J 


This represents the sum of heat energy H and kinetic 
energy ox fluid per unit mass, A knowledge of the time rate 
of change of the stagnation enthalpy may be gottam by the 
following procedure® 

From equation Q I ■ / have 

„ 


Pr 

J)hv 

~jn~ 


>L Jl C 
-I? + Mr 


_J>H 


-+ i 


pJt~ 



be written as 


i®e e dQ — © and we have 


and (. >' ,y Jare the energy equation for the 
il adiabatic flow of comprassible flow® 
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EQUATION OF state » 

The measurable quantities of a compressible 
substance are its pressure p, density f , and temperature 

T * Xt is found that these quantities are connected through 
a functional relation of the £orm<> 

p (P# / , T) =0 

Where P in a single valued function of the 
variables p, I , T such an equation is known as the equation 
of state of the substance* The form assumed by F depends 
on the nature of the substance® 

.that the pressure in a fluid is ascribable to 
the momentum change produced in its molecules on impinging 
against a small plane rigid surface. The density is prop- 
prtxonal to the number of molecular structure of matter. 

In addition, since the randomly moving molecules of a 
substance have a mean kinetic energy associated with them, 
this leads to the notion of, temperature, being defined as 
a measure of this kinetic energy of random molecular motion. 

For the moment we confine our attention to 
gaseous substances. In some cases, the molecules of a gas 
have but negligible volume and there are vi rtually no mutual 
attractions between the individual molecules. Such a gas 
is said to be a perfect gas and its equation of state 
assume the simple form. 

P = R f T 

Where R is a constant for the particular gas 


I 

J Vv. 

1 




under consideration 


I __l ; 
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be some thermodynamic is function of a substano 
expressible in terms of any two of the measurable quantities 

have seen that this is always possible whenever 
the equation of state of the substance is known, it may 
happaii tiiac. x X we express as a functioning, say, p and 

V ( V 1// f ) 9 then in changing from an initial state A to 
another state B the value of depends only on the condi- 

at A and B and not at all on the paths joining them* 
in these circumstances, ^ is termed a function of state. 

Suppose that the differential of f is expressible 
in terms of d R and dv by the relation. 


If is independent of the path taken in 

changing from the conditions at A to those at B, then 

must be an exact differential. The necessary and sufficient 

condition for this is 


^ )P 




V“ 


then <=fc (p, v) is a function of state 


The first law of thermodynamics gives 
dQ = p dv + du 

for a perfect gas, du ~ T 

q[U = (_ ^/p-J cl ( ^ ) 

l=>^ ~ 97 


using the equation 
Thus for such a gas 


/ 


identity dQ with d<4> , we have 

N cw 


( 



J h 





<si nets Cjo Cv # dQ is not an exact differential. This 
means that the quantity of heat Q added to a unit mass of 

fJU3 not d function of state, it depends on the both 
chosen in the change from a given state A to a second state 


Using pv - RT the relation for dQ for a perfect gas 
may be written in the alternative form 


oK = (^Jd- 't (±-}<Up 

where cl S — d 8 j~f 



so that 




s called the entropy per unit mass 


Flour for which S is constant! 


and 


Ir a change is made so that t h€ ontftpy of every single.' 

•V 

particle of the working substance stays constant, then such 
a change is termed isentropic when the entropy of every single! 
quantity of a substance of fixed mass is the same and stays 
constant in any change, the change is said to be homentropic. ! 




“‘“ UCtin9 “-Possible fluid could 

>mit s ound waves which ard 1 

Ch dre lo ^tudenal in kind, no 

’ ran “' SrSe pr ^gatica in possible for such a 

• More over an inconpr. sslb i a fluid. „ h i ch a non 

=ting can not sup pert aoy kind Qf uav£ „,. tion #t #u 

flUid iS however transverse wave 

1 thr °ugfa it is possible when ever * 

Ae v ' nen ever a magnetic field 

ssent Til [ o 1 


M ^ multiply both side 
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sagjgg mm, 

mm- 


j £j 




= ( ur j 

c^ 


L<atI L Uly^I yj 
where 


H <crjj 


through vector identity 
~X7 ^-jJ( ■=_ ^orcus| <sLv H 

V 2 H - — Cu^j Ct-irJ W* 


(Lm'tJ (_MrJ 


— *) — 


<W (vUH) + l "j ^» 2 H 


We can written another way 


C^rKv*^ 


«J) \ ~' 5 > 

V* (_ K ^ / V" 

— ( xA ^7 .) H 


- ( ^ ^ j fr*- 


if the fluid is incompressible ^Ijy \/ =- o 

■ CurJi — f & ^ ^ v ; h’- 


C f? w v* 


<?h 








This gives the equation of magnetic field for incompressible 

coni' rti nc fluid 







viscrus 



snocK to be a surface of discontinuity 
velocity, pressure, density and temperature of the fluid, 

shall find the relation between these quantities in front of 
and behind the shock* We choose the coordinate system, such 
that the shock is stationary and the fluid moves through it. 
Under these conditions in flow is steady we assume that the 
fluid is an ideal gas so that the perfect gas law holds both 
m front of and behind the shock. ■ 

SHOCK 


U 2, P 2'X 


To find the velocity u 2 , pressure and density behind 
the shock, if the corresponding values in front of the shock 

are given. 


Equation of continuity 


Equation of motion 




This is known as the Rankine 


.now that ratio of 


ave from isentropic condition 






Lu — -U. c= Cv / ? '(y'-fO / ^ 


1 2 e spee® of sound 

shock wave just after the shock wave and 
velocity of sound { m . J., ' q « a 

h* __ C\jjc ^ ( y-t-) J 
-f, — ~ 


Substitute the value of (1*18) and 
ion (1.17) . 






Wave J U l' u 2 Wl11 never increase to the critics 
sound* Let and Mj denote the maeh number j 
after the shock wave. 


Since by the R anki ne-Hugoni te relation. 

K _ Cv->H " t r+,) 
b, C y-' 1 ) fz ~~ ( 


r-> 17 

(_ -T+ K 
C -r- 


from conservation law 


Substitute it 


Uz. 

~~ . Vv\ 

Lu 

f iii- 

y~\~) 

. ^ 


if A- 







divided it 
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GENERLIZED RAftKINE-HUGONIOT RELATIOt 

" shock WAVE 


1EGNE 1 ] 


Shock wave under the effect of finite viscosity and 
thermal conductivity when the gas is conducting,, When we 
consider a shock wave of a conducting fluid the surface of 
discontinuity is very thin and for a non viscous having to® 
thermal conductivity* 

In the case of conducting shock wave under the effect 
of finit viscosity the surface of discontinuity has certain 
thickness know as shock zone* We consider a stationery shock 


zone of thickness in tube of uniform cross section. The 


flow is one dimensional and shock zone is separate in region 


1 and 2» 


Changes continuously from -P. , p, u 1 to P_ 4 u~ Thi; 

• 11 2 ^ ^ 2 ® 


continious change from the shock zone is give to finit 


viscosity and thermas conductivity of gas* 


Now we consider the motion of element of fixed mass 


moving along with the fluid particle with velocity u 


the flow is one dimensional, 


That there is a magnetic field of intensity H = H(x) j 


at right angle to its motion where H varies continiously 


from H. in region 1 to H , in region 


Then the pressure of magnetic field will increase the 
pressure by ( The magnetic energy per unit volume 

of the medium is ( <1 ) o 











I 

1 ■ ■ 
B . 


| 

ft 

t 


r C|li 








J 3 ' , iFf ; Lj, Z OAj£T I 


conservation of mass remains unaltered 


The momentum equation 

Pi *+ -p, 1 


xs modified to 


oj -?i vj ( 2 t- ( ^ "* ( 


we obtain the energy equation 


ii, f V*-U. ” 

1 Ctf>1 
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To account for the effect of the magnetic field we 


must modify the effective pressure from p — v' 


£ ^ jdnj 


S v. ' S 

to ( ^ Ari2 'l<3K J ^ ^ and also to the right 

hand side we must add the rate of gain of magnetic energy of 


the el erne 


nt. Now the magnetic energy of the element of mass 




Hence the rate of increase of this magnetic energy is 

4f (Af i*.) I )? - ("ixT J 

1 . . e . / Un 2- \ 

CO. /ffM dr*- ^ y 

Hence the energy equation becomes 

4^ =■ ^ <& 0 M + ^ p- 

LA yU- ° 

4* ^ ^ * ca c 4 /_£ ) -t ( f ^ <W L * 2 ' 

Cr-i J <**■ < ' 


/) 4-<f x Sn + 




l T- 




JU H 


Jividing through by f u A J* . gives, since “f M == constant 


— ? — -A + u. J_ u i - ia sil - — d Z 

a a j | z 


integrating through the shock zone and using du/dx = 0, 


dT/dX = 0 in region 1, 2 


V _£l. 

fr h 


Jj H> 

L-i 'K f j 


~f_ 1^ 2. 

cr-V ^ 
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We thus have three equations for the four ratios 


Thus 


Integrating this 
through the shock zone. 



Hence that one of the effects of the magnetic field 
is to decrease the strength ( P 2 - P 1 ) / P 1 of the shock 
also, for a non-conducting gas the velocity in region 
must be greater than the speed of sound a, of the media but 


with the transverse field present for a donducting gas 


Shock formation requi red 


momentum and energy and maxwell 1 s electromeganetic equations* 
Several types of discontinuities can exist in ideal electrically 
conducting fluids in the presence of maganetic fields. The 
discontinuities characterized by the condition that both the 
mass flow and density change across them are different from 





SIMILARITY CONSIDERATION! 


The motion of a gas or a liquid is said to be one- 
dimensional, the only possible one dimensional motion are 
produced by spherical, cylindrical, and plane wave® The 
simple wave solution is limited to plane wave moving in to 
a non uniform region® However, the problem of plane cyline- 

drical or spherical. Symmetrical shock wave moving in to a 

% 

non uniform region are more complecated. A fairly general 
approximate theory for shock waves can be obtained, known as 
similarly solution® [ 7 ] 


men tzne characters tie parap^er include two constants 
with independent dimension in addition to r and t, the partial 
differential equation satisfied the velocity, density and 
pressure in one dimensional unsteady motion- of a compressible 
fluid® Can be replaced by ordinary differential equation for 




The solution of these ordinary differential equation 


can some times be obtained exactly in closed form. Such mot 


non heat conducting, so that the motion does not involve any 
kind of physical or chemical change 


motion, continuity and energy take 


The equation of 
the form 


OJtr 



~t 


i £ b 

T 


D 


(1.27) 


it 

o )Jr 


+ 


$ u 


-+ 


[ y-t ) 


-f VA 


n c 


(1.28) 


It?) ■+ J=-° ci.29) 

/ 

Where is the adiabatic exponent? Vsc o for plane 
flow, Y ) f 01 flow with cylindrical summetry ; and. J?. 

for flow with spherical synimetry. 

which all the flow quantities take the form t m f (r/t n ) . 
These have simplifying feature that the partial differential 
equation reduce to ordinary differential equation with indeper 
dent variable r/t n . These similarty conditions originated, 
from dimensional analysis is based on the fact that the only 
parameters in the problem are E with dimension ML^/T^ and 

with dimension M/L the only parameter involving dimension 
of length and time is with dimension LS/T^ or some 

function of it. Taking this is to consideration we can drive 
the various quantities that arise in solution f 7 I 



THERMAL RADIATION 


Radiation phenoraana have acquired interest for gas 
mairly since attention has been attracted in Science and 
Technology to such phenomena as nuclear explosions, hyper- 
sonic motion of bodies in the atmosphere, powerful electric 
discharge and astrophysics! problems, the modern trends of 
aerodynamics are towards high speed and high temperature 


as well as low density and high attitude* At higher re-entry 


speeds such as that for a mars probe, the interaction of 


aerodynamics and radiation fields becomes important radiation 


gasdynamics is conceived with the study of the effects of 


thermal radiation in a very high temperature gas flow. In 
radiation gas dynamics, the differential approximation for 


:ay gas of arbitrary opacity is largely used* 


The theory of radiative transfer and radiant heat 


exchange was created and developed in order to understand 


processes which take place in stellar media and to explain 


the observed luminosity of stars. To as large extent this 


theory can be also applied to other high temperature system 


considered in modern physics. The thermal radiation is 


characterized by frequency ^ of oscillation of an electro- 
meganetic field or by the wave length fa related to the 
frequency and the speed of light c by the relation ^/"P 


If a fluid body is heated non-uniformely or if energy 
is released within the body, a thermal flux transported by 
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diffusion and temperature equilization, for a weak shock, 
when heat conduction is present but viscosity is absent, 
the gas can not make a continuous transfer from initial to 
final state. A discontinuity is necessarily formed, which 
corresponds .to a viscous compression, shock® 

If the heat conduction flux is proportional to temperature 
gradient than all flow variable with exception of temperature 
undergo a discontinious jump at the discontinuity and an 
-isothermal shock" which are caused by electron and radiation 

heat conduction. 

However, for extremely strong shock waves, when the 
energy, density and the radiation pressure become sufficiently 
large in comparison with the every and pressure of the fluid® 
The situation cfoa.fg@z and the ges in the... shock wave makes a 
continious transition from the initial to final state through 
radiation conduction alone [ a] (sj [ 10 1 

(i) radiation pressure 


Which may be expressed as 


■R 


1 

3 


T ' 


( 1 . 30 ) 


Where T is the temperature of the gas and s al is known 
as the Stef an-B o 1 tzmann constant® Total pressure at each 
point of flow field is equal to the sum of gas pressure and 

radiation pressure. 
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(il) RADIATION ENERGY DENSITY 

The radiation energy density per unit mass of ' the 
fluid is given by 


tzp. 

~T 


< 5 T 


f 


Cl. 31) 


Where f is the density of the fluid. 


(iii) RADIATION FLUX F 


The radiation flux F is given by the formula, 


P 


ac 7- 

~ ~T~ 


(1.32) 


Where c is the velocity of light in vaccura. By taking 
into account the effects of radiation and considering the gas 
to be inviscid and non heat conducting. 

The equation of motion and energy 

i - - ‘Irad ( ftffi.) (1.33) 

DJc Q ■ 


and 


f . 


cl* 


,v (H) oL ' jf 


where 


-e 




D , J_ a * 


(1.34) 


(1.35) 

and is the total energy per unit mass 



of the luminescence of a tints ocuj. 
•x circle is boundary of the lumino 
miter circle is the shock front® 
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VISCOSITY 

Consider in the real fluid* the surface of discontinuity 
should be replaced by transition region. If the thickness of 
the transition region is very thin so that we are interested 
in the flow inside the transition region or in flow phenomena 
closely related to the transition region? we must consider the 
effect of viscosity. 

In many problem of fluid flow, the fluid is actually a 
mixture of several gases, i „ e. air is a mixture of oxygen, 
nitrogen, and other gases. If the concentration of the gases 
in the mixture remain unchanged during the flow field, we may 
consider the mixture as a single fluid. 

Viscosity represents that property of an actual fluid 
which exhibits a certain resistance to alteration of form. 
Although this resistance is comparatively small for many practi 
cally important fluids, such as water or gases, it is not 
negligible, for other fluids such as oil, glycerine etc. , this 
resistance is quit large. In a viscous fluid, both tangential 
and nortttil' Zonae of -the HiHeiic energy of flow 

will be dissipated as heat through the viscous forces. 

First introduced the concept of the boundary layer by 
stating that the viscous effects are confined to a very thin 
layer near the boundary. The viscous effects are confined to 
very thin boundary layer whose thickness tends to be zero as 
the coefficient of kinematic viscosity goes to zero. The 
Navier-stoks equation can be simplified by such consideration. 



analysis so that equation of motion can be transformed to 


numer. 


suggested by Richtmyer and Von-Nenumann is given as 



- ^ } 

d. 

J 

dX) 
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SPHERICAL SHOCK WAVES IN VISCOUS MAGNETO GAS 

. "dynamics 


***** isl&'ik 


The resistance of shock waves gasdynamics flow field 
introduce free boundary discontinuities in to Physical Parameter 
of the system. A method for avoding such difficulties, particu- 
larly for the numerical calculations, has been developed by 
Richtmyer and Von Neumann. 

They observed that the addition of a particular viscosity 
like term into gas dynamics equation could lead to continious 
shock flow field in which the finite thickness of discontinuities 
at theshock wave was removed and replaced by a region in which 
the Physical Parameter change rapidly, but smoothly. 

In the present work, the artificial mechanism of viscosity 
in the presence of magnetic field to smear out discontinuities 
of the physical parameters from the flow field. In order to give 
a meaning to the other wise physically unrealisable magnetic 
field with the spherical symmetry, the magnetic field, is replaced 
by an idealized field such that lines of forces lie on a hemi- 
sphere whose centre is the point of explosion. We have used the 
Runga-Kutta method to obtain numerical solutions in viscous and 
inviscid regions. We have shown that field variable change 
rapidly when the magnetic fiS^d is impo'sed in both the viscous 



-A. 

and the non— viscous regions© 

The equations of motion of a fluid having infinite 
electrical conductivity with artificial viscosity and expressed 
in spherically symmetric form, are, 

JF ^ u fk *+ f i 


JL 

o t 


-+ M 


d A~ 




< 5 ^ 


-f- 61 ^tL 
&Jr dX 


-f H 


J M 


d* 


A* H 


-■"D 


( 2 . 2 ) 


(2.3) 


and. 




-f u 


<#> 


^ fj? -h i r-ljtj 


f 


+njL 


(2.4) 


Where u is fluid v exocity, f be density, h be component of 
transverse magnetic field and q is artificial viscosity. 

We seek the solution of equation in the form 


p (%,+) --- p-fyx I ( 
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and the non-vi scous regions® 

The equations of motion of a fluid having infinite 
electrical conductivity with artificial viscosity and expressed 
in spherically symmetric form, are. 


it. -f- u -+ 4- 
df 1 


? M ^ 


^ u 


}Jr 


d h- 




a. c 

a-A- 


+ Mb 

U(2.l) 


0 

/ ^ 

-h 

3-H' \ 
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J — 10 

(2,2) 

-t 

, H 

dX 

- t 

>mh 

~"'" D 

(2.3) 


and. 


?jp 

^ 4 


■f U 


it. 

^ /}?•+•( J 

u + 

1 

<) 


i 

( d/k 

f 


( 2 . 4 ) 


Where u is fluid velocity, ~f be density, H be component of 
transverse magnetic field and q is artificial viscosity® 

We seek the. solution of equation in the form 


jo (h, +) — U P 


j fl %,r) 


A 


fa 

£> , r 
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and 


Where x 


R 


j/r 

i± i 

d-t' 
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( a., p) 

- X 

i-o ^ _i2 
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n, j 


* R 

being a function 

of time only 

we 

assumed that. 
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Di £ £ ®it®b tistion with ir^op©ct to jr 
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( <H ) C I'* (3 ) ( \-) Cl) represent 

the equation in the viscous region 'X- ~?/ | 
for the region without viscosity £ © £ x < \J q = 0 i*e 
G = 0 the equation reduces to* 
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The equation are solved numerically for the viscous region £ X >t ) 


The jump conditions are given by 
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CONCLUSION 

The equation (H) , (I) , C J) , (K) and (L) are solved 

numerically using boundary conditions for the viscous region 
( > ) for which f 7 ’ «c p and the numerical solution of 

equation (M) , (N) , (0) , and (R) are obtained in the region 

without viscosity ( o < X < I ) 

We conclude that the field parameters change rapidly 
in the inviscid region as well as in viscous region when the 
magnetic field is imposed. The variation of velocity, pressure, 
density, magnetic field and viscosity for both the viscious and 
non-viscous region has been illustrated through the graph 1 to 9 . 
We also infer that in viscous region, the field parameters 
attain a maximum value at the shock front when the artificial 
mechanism of viscosity is weaker; i.e. for smaller value of 
K *= 0.0349 and starts decreasing rapidly was we move away from 
the shock front except the pressure which decreases rapidly 
only for a narrow region and the increases instantaneously. 

The magnetic field has the significant effect on the flow 
parameters for this value of K but for K » 0. 349. 

The effect of magnetic field does not illustrate any 
significant change in flow parameters. This shows that the 
effect of magnetic field does not play any important role when 
the artificffiial mechanism of viscosity is stronger. 
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SELF-SIMILAR CYLI NDRICAL SHOCK WAVES WITH RADIATION HEAT FLUX 
INTRODUCTION 

The problem of propagation of shock waves in a non— 
homogeneous medium is of great interest in exploring the 
effect of explosion in the stars and atmosphere of the earth. 

The solution for cylindricaily symmetric flow has been 
obtained numerically by Lin (19), Ray (20) has discussed the 
problem of point and line explosion and found an exact analytic 
solution. Analytic solutions in the three cases of plane, 
cylindrical and spherical flow have been noted by Sakuri (22) , 
Rogers (21) has also studied the similarity solution for these 
three cases in uniform atmosphere. 

In the present paper the problem of explosion along a 
line in a gas cloud has been discussed. Similarity solutions 

h<3v4; b&oK descnbtv^ bM. prop^&tfoVt oj- Q cy yuhic&Jj 

shock in non- uniform atmosphere taking counter gas pressure 
and radiation heat flux into account. The radiation pressure 
and radiation energy have been ignored. The gas in the 
undisturbed field is assumed to be at rest® We also have 
assumed the gas to be grey' and opaque and the shock to be 
transparent and isothermal. 





Put the value from the equation (3.4) , (3.5) and (3.6) in 
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For an ideal gas 

P = 5"! T 
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Taking Rosseland* s diffusion approximation we have 
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isothermal shock and the conditions are 


£ - e, 


^ 5 U 


? rr ^ La i_ "t* 


P, v 
v p, 


r rh 


2 — Vvx < 




vi f ' 
V >>1 


Substitute the value of Vm ^ | (v~U^J - j(? \j „ 


■»- [ v/ " L '0 -y' f'l' 2 - f, v/ 


V 


r-Jv 


y ho' 


























CONCLUSION 


The differential equations (a) to (D) are numerically 
solved by Runge-Kutta method. The numerical results for a 
certain choice of parameters are reproduced in the form of 
table we conclude that the field parameters decreases in the 
radiative shock front® The radiative velocity (u) * density 
( *f 1 * Pressure (p) and heat flux (q) all are decreas in 
first and second set but shock front is increases in table 
first and second. 

The table of certain choice of parameters with reproduce 

in the form of tables. 


First *ct 




M 2 * 20, .V * 10, w = 

- 1.5, * = | 


n 

u/u 2 

• p/fe 

Pil>2 

«/<?; 

1.00000 * 

1.0000000 

1.0000000 

1.0000000 

1.0000000 

1 .11 10 (JO 

0.9735775 

0.8409035 

0.7628222 

1.045415 1 

1 .200 00 

0.95301! 1 

0.7323471 

0.595 1056 

1.0873912 

1.30000 

0.9369498 

0.6567086 

0.4703231 

1.1243515 

1.40000 

0.924 2489 

0.6043024 

0.3731744 

!. 155571 0 

1.50000 

0.9 1 3 8.su 6 

0.5706075 

0.2944013 

1.1806336 

1 .000 00 

0,9044285 

0.5 u» 42 2 9 

0.228224 2 

U 99 056 2 

1 ‘OUt hi 

0 K9 U»39H 

0.5749437 

U, 17 12666 

1.209619 8 

l .81)000 

0,874 1 M 1 

0, 706400 36 

0.1233300 

1.210213 1 

i. you oo 

0.8198924 

2.9936543 

0,0742085 

1.1910570 

1.01000 

0.8036516 

33.050653 

0.0413089 

1.184851 5 


TAHI.li 11 

Second set 


y-lM 2 

« 20, A' - 100, hi - 

- 1,75, a “ 7 


n 

W«2 

P/P; . 

PlPi 

0/4; 

1.00000 ■ 

, 1 .000 0000 

1 ,0000000 

1,0000000 

1,000 000 Q 

uOoou 

0.96! 1185 

0.8255005 

0.8181402 

0.9282533 

i. 200 00 

(1.928 0H-1 2 

0.7043508 

0.6924694 

0.8824891 

1.30000 

0.899762 1 

0.6169417 

0,6020795 

0.8529418 

1.40000 

0.8752778 

0.5518644 

0.5349040 

0.8339617 

1.500 00 

0.853 9-, u 3 

0.5021333 

0.4836044 

0.822 1112 

1.60000 

0.8352256 

0.4632864 

0.4435191 

0.8152099 

1 .700 00 

0.8 18 6828. 

0.4323696 

0.4115770 

0.8118257 

1,800 00 

0.8039696 

'0.4073669 

0.3856908 

. 0.8109901 

1.90000 

0.7908040 

0.3868654 

0.3644028 

0.8120310 

2.00000 

0.7789563 

0.3698514 

0.3466707 ' 

0.8144719 

2.10000 

0.7682379 

0.3555831 

0.3317335 

0.817968 2 

2.20000 

0.7584934 

0.3435072 

0.3190253 

0.8222664 . 

2.30000 

0.7495935 

0.3332043 

0.3081179 

0.8271774 

2.40000 

0.7414304 

0.3243518 

0.2986825 

0.8325582 

2.500 00 

0.7339134 

0,3166982 

0.2904632 

0.8382994 . 

2.60000 

. 0.7269657 

0.3100448 

0.2832582 

0,8443164 
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Abstract. The point-source, spherical magnetogasdynamics shock wave moving into a constant density 
y-latv gas is considered in the limit of infinite shock strength, from the point of view of the Richtmyer-Von 
Neumann viscosity technique. Numerical solutions of this problem has been obtained in viscous and 
noil-viscous regions. A similarity solution of this problem is shown to exist. We have shown that field 
variables change rapidly when the magnetic field is imposed in both the viscous and the non- viscous regions. 






1. Introduction 

The existence of shock waves in gasdynamic flow field introduces free-boundary 
discontinuities into physical parameter of the system. Such discontinuities cause 
considerable analytic as well us numerical complications in the treatment of gasdynamic 
problems. A method for avoiding such difficulties, particularly for the numerical calcu- 
lations, has been developed by Richtmyer and Von Neumann (1950). They observed 
that the addition of a particular viscosity like term into gasdynamic equations could lead 
to continuous shock flow field in which the finite thickness of discontinuities at the 
shock wave was removed and replaced by a region in which the physical parameters 
changed rapidly, but smoothly. Following this method, Lax (1954), Latter (1955), Rrode 
(1954), Colgate and Johnson (1960), and Christy (1964) solved various shock problems. 
Sachdev and Prn:.ad (1966) us -d the method of artificial heat conduction to sm~ar oat 
the shock in ordinary gasdynamies. 

In the present work, we use the artificial mechanism of viscosity in the presence of 
magnetic field to smear out discontinuities of the physical parameters from the flow field. 
In order to give a meaning to the otherwise physically unrealisable magnetic field with 
the spherical symmetry, the magnetic field is replaced by an idealized field such that lines 
offerees lie on a hemisphere whose centre is the point of explosion (cf. Summers, 1975). 
Wc have used the Runge-Kutta method to obtain numerical solutions- in viscous and 
inviscid regions. The numerical calculations have been done on a DEC-system 1090 
computer installed at l.T.T. Kanpur by RKGS programma. 

2. Formulation of the Problem 

l The equations of motion of a fluid having infinite electrical conductivity when expressed 
■in spherically-symmetric Eulerian form with artificial viscosity term as suggested by 
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The existence of shock waves in gasdynamic flow field introduces free-boundary 
“ ; ; A method for avoid, * such divides, pardeulari, 
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Richtmyer 


field, and artificial 
taken to be unity 


atien of state is assumed 


where 


The form of q chosen in (7) has been made from the point of view of admitting similarity 
solutions, since all forms of q may not fulfil this condition. 

We assume the existence of a similarity solution of the form 

p(r, t) = p 0 R~ p f{x) , p(r, i) = p 0 >K*) - (8) 

fi(r, t) = R~ x <p(x) , H{r, t) = »;(*) - 

where ,v = r/R, R being a function of time only, p 0 is a 
dimension of density, and a, p are constants. To simplify 


m orbiirary constant having 
? the subsequent calculations, 


SHOCK. WAVES IN MAGNETOGASBYNAMICS 
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it is assumed that a = -*/?. We change the independent variables from (r, l) to (.v. 0 with 
the aid of the following relations and later pul £ * t\ we find that 

d d xR 1 d ' 

it ~ it R 8x ’ dr R dx ’ 


where 


Equations (1) to (4) then take the forms 

f + §' + tin' n 2 _ a 

-A(otcp + x(p') + (p(j> + “ + ^ 


-Ax) + W+~” °> 


ip'(<p 


rj' (<p - Ax) + rjltp' - a4 +~J - 0, 


4> 




where a prime denotes tie dSre-ir.ioh with respect to * and the quantity 8 * «*“«* 
to q by the relation 

gix) (13) 

R 


l^X 2 f - ^ = ' r,2a 


It is assumed that 

R a R = A (a constant) . 

If E be the total energy, then 


“ 2 ' + 1 dr , 


r(p T + y-l 2 


which becomes in the terms of present similarity variable as 
E = 4 ti A 3 2a rt> 

Since R is a function of time only, the energy £ will be independent of time if « - t- 
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We now introduce new dimensionless quantities 

f r 8 P = t N = - • 

F = - , G = r a A 

A A A 

. ^ „ /on tr, tin and (15) with the aid of (16), we get 

Transforming Equations (9) to (13) ana 


( 16 ) 


and 


F' + G' + Ntf N 2 _ o 
-(ip'+xn + pp'* — + j x * 

(17) ’ 

% 

p 

\p'(P - .r) + $P' + 2 iA — = 0 , 

(18) 

It# 

,V</> - x) + ;v(p’ + LD’ 0, 

(19) 

iti'd 1 

_, F ^P- x) r- y FtSyfF G (p-xW-o. 

.// 

(20) 


g = — 2 .v 2 r (|F’| -P r ), 

2 

(2!) 

TlL- 


(22) 

1L 4 


Tavlor (1950) has shown the existence of a solution of equation of the type 0 7) to (21) 
f or a diverging flow from a point source with a shock discontinuity at -v - 1, since fo, 

such a diveiging flow P > 0 at d Equation (21) will he divided separately 

The treatment of Equations (1 /) to (2U) ana ' * 

for the region where (he viscoshy is absent and the region where it ,s present. First the 

region without viscosity will he treated and. thereafter, the region ,th viscosity will be 
explored. They would later be matched through numerical integration. 

3. Solution of the Equations (17)-(21) 

3.1. In viscid flow 

From the boundary condition for a diverging flow field it follows that, at r . R, the 
gradient of velocity is positive and, therefore, the viscosity term ,s aero - i.c., C - 0. The 
equations defining the flow and the field equation reduce to 


F' + A IN' N 2 „ 

.(I p + .,n + w+— -v 


(23) 





SHOCK WAVES IN 


MAGNETOGASDYNAMiCS 


1 

( p - x)ifr' + + 2ij/ — = 0 5 




We use the boundary 


X V A V ' 

- 3 F+(P~ x)F' - jH p ~ x ) ^ 

Rendition to obta* the numerical so.ution of the 


,.2. Viscous H.ow it i s essential to have discontinuities in the distribution 

In order to introduce viscosity, it nosition of these discontinuities at r - R, 

of flow variables. We may arbfirar V the visC ous region *r 

so that the problem now remains to obtain 

which ^ 0- 

Equation (21) may be written as 
1 


KX . ft and so by Equation (27) F = 0- This condition 

at the transition point x - 1, G > y ^ transit ion from x = 1 onwards, 

determines the magnitude of jump in the slope 


4. Results and Discussion 


To illustrate the behaviour of the present 

Equations (23)-(26) are oblame , numerically for the viscous region (sc 2 D- 

r— 


/•(!)= U~ 6 )’ /, ' (1) = (l " d) ' 


N( 1) = 


• u a on/t inst behind the shock front, i.e., Po/P* d 

where d is the ratio of densities just a ® a The ralio of densities - i.e., pjp u for 

and Alfven’s Mach number - vW ■ _ „ 0349 lt » 4.054 and for K - «■» 
v - «*<£ Related our problem for these values m the 

presence of the magnetic field. . _ n (hc invisci d region as well 

We conclude that the field P® raine ' ‘ )cld is ilu posed. The variation ot velocity, 
as i'n viscous region when the l »« bo th the viscous and non-VlSCOUS 

pressure, density, magnetic held, dad vtscos J ^ g Wu ^ infer that in viscous 


nmtm 
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Ato'i. sumo ll “ *" a 

jsl. - »«« - “ p "“ e 

falls us the inverse square ol die distance 


1. Introduction 

The problem of propagation of shock w "“ “ Mmosphere of the Earth, 

interest in exploring the effect o exp tost obtained numerically by Lin 

The solution for eylindncall, aymmetn Sow ‘ fa explosion gnd found 

cylindrical shock in non-uni orm r e and radiation energy have been ignore ■ 

heat flux into account. rad >“ P j b at rest . W e also have assumed the gas 

of the explosion is constant. 



2. Self-Similar Formulation 




behind the wave are 


^3-268. 0OM-640X/84/ 1022-0263$ 00.90. 
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conditions, following Singh and Vishwaknrmn (19R3), arc compatible when 
km co, A = 0 , a = - (4 + co) , b = 2. 


n - - 2 , P 


- 1 and a - 


Cl) + 1 


69 


U> 




rv 


5- 


Hence, the pressure distribution becomes 
/h = BR~ 2 ■ 

F/!-;at : ons (1 }— (3 > and (:) ar "< ! 'ii o • 0 

.fl'{rj)IA _ >?(4 + (o)v'(*l) ~ (2 + to)v( n) ? 
0(5/1 f/ [2-(4 + coM^)] 


j" | ,-VSV 


r(f?)M = nO?)A4 

2 + o> 


{2 - (4 + co)i>0?)} + Kr?) (i- 0;)_2J )1 
t]( 4 + co) -1 


+ 


t]( 4 + co) 


f'Oj)/a _ P’(n)!A 


F{r\)jA F(>i)IA 

. P(n)IA 


2 - ( 4 + (,J M^ 
L (>’ - 1)(4 + cj) 


] 




(20) 


( 21 ) 


( 22 ) 


(23) 


f(co + 2y + 2W») - y(4 + + } 

rj(y - 1) (4 + co) J rj 

f p'(n)M _ 


(24) 


FWM w (*W /2 fo _ ( 4 + co))? 1 6-^- - 

^■-''SSS^L 2 * n WM JJ.W) 


(Q(f?)M) <3/2 


I 


when 


Nm 


4<tCju 0 
3 ¥ 3/2 A'- a 


= a non-dimensional radiation parameter , 


•y* 

I (F(>mm>m (V2 ^ [2 _ (4 + w)v{>1) ] _ 2 

N {P(n)IAV 2 


Pin)h 4 


X 


(26) 



rc/i + <•.,)„( n ) + - [2 - (4 + 

r ~ x2 P( r?)/xi 1 
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* V 


conditions, following Singh jind Vishwakarma (1983), arc compatible when 
k = co, A = 0 , a = - (4 + co) , 6 = 2, 


n = - 2 , fi ; 


• | and a = 


w + I 


CO 


IV 


to" 




Hence, the pressure distribution becomes 
P\ = BR~ 2 ■ 

Ivv.at'ons f!)-(3 ; } and {:) «•: no Inu I to .«*. f > -r* 

p.'{rj)/A _ tj (4 + co)t/(r?) - (2 + co)r(Q ^ 

"'nOti/A ~ tj[2-(4 +«)»('?)] 

P'(n)IA _ Q(>l)IA f r)v'(tf) {2 - (4 + co)o(rp} + vjrj) 0 : 0l) l_0 
’^j/7 ~ P{n)l A L »j(4 + ©) 

2 + W 


r/(4 + co) 


F\n)!A = r 2 - (4 + a) » (>>)' 

~F(n)IA F{n)IA (y - i)(4 + o>j . 


( 2 ( 1 ) 


(21). 


( 22 ) 


(23) 


P(q)/A r (co + 2y + 2M tj) - y(4 + co)r;r'(r;) ~ 2 




F(t?)M L 


n(y - 0(4 + co) 


+ - , (24) 

n 


when 


.mi a = _ N 


P(rj)IA (Q(tj)/A) (3/2 


- *) 


2 - (4 + co)f? 


p’(tj)/A _ »' (>om V 
p(n)/A n(>j)/A }. 


Afm 


4 (rCfA() 


= a non-dimensional radiation parameter , 


, (25) 


(2b) 


3r 3/ 2y4 l-a 

1 <KW) mmr" « _ (4 1 

n {P{nW 2 m)/A 


;■ I 


rot + o>)nf »;) + re) - [? - (4 + coMrpI 00?) (v(tj) - 0] 

' d 


* Pin )/ a 


{2 - (4 + co)o(r/)} 2 - (4 + co) 2 ----- 
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The approximate shock conditions are 

i r* t 1 
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FOl o) _ I f 

^4 2V4 + «y 


(28) 

(29) 

( 30 ) 

, (3D 


_vr_i — iT. - ' " " 

-"T - - oU + coj ly 2 M 4 J 

• i V ilues for our numerical calculation, where we assume that 

which are the imual values 10 


3 . Results 

toro “ 

v i/C4+«>)jp(r)) 




and 


U 2 V 112 

u'VW/ 

/.. \u)/h + “') a(ii)/4 

M® 

a(t? o)M 




w) /(4 + w) P(i])jA ' 

p = (2s\ 
'p2 W 

p(iio)M 

/ \i3 + »V(4+“>> F(f?)/-4 

<? = 32 
fl, w 

FOl Jo)/- 4 


i' 1 


(32) 

(33) 

(34) 

(35) 


^ .Nature of the field variables mayo lers; 

tables. Nature following two sets m v 

calculate our results for tn = _ L5 , 

, jr9 . OA ?V = iU , 


7=3 


'■ u * = 20 , * = 10 , 

M 2 = 20 , # = 100 » 


to- -1-75, 


1 

(X * 3 * 
3 

*<* = 7 ’ 


w ]V = 100 5 tw - 

sssSSssssskw- 
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TABLE 1 
First set 


y = M 2 -20, N = 10, « 


- 1.5, * = 5 


t. 000 00 
1.10000 
1.200 00 
1.30000 
1.40000 
1.50000 
1.60000 
1.70000 
1.80000 
1.900 00 
1.91000 


1.0000000 

0.973 577 5 

0.953011 1 

0.9369498 

0.9242489 

0.9138366 

0.9044285 

0.8936298 

0.8741511 

0.8198924 

0.8036516 


1.0000000 
0.8409035 
0.7323471 
0.6567086 
0.6043024 
0.5706075 
0.5564229 
0.5749437 
0.706 46636 
2.9936543 
33.050653 


1.0000000 

0.7628222 

0.595 1056 

0.4703231 

0.373 1744 

0.294401 3 

0.2282242 

0.1712666 

0.1233300 

0.0742085 

0.0413089 


1.0000000 

1.045415 1 

1.087391 2 

1.124 351 5 

1.1555710 

1.1806336 

1.1990562 

1.2098198 

1.2102131 

1.1910570 

1.1848515 


TABLE II 
Second set 

20, JV- 100, w = 


1.00000 
1.10000 
1.20000 
1.30000 
1.40000 
1.50000 
’ 600 "O 
I 700 UU 
1.80000 
1.90000 
2.00000 
2.10000 
2.20000 
2.30000 
2.40000 
2.50000 
2.60000 


1.0000000 
p.< 61 1185 
0.9280842 
0.899 762 1 
0.8752778 
0.8539463 
0.8352256 
0.8186828 
0.8039696 
0.7908040 
0.7789563 
0.7682379 
0.7584934 
0.7495935 
0.7414304 
. 0.7339134 
0.7269657 


1.0000000 

0.8255005 

0.7043508 

0.6169417 

0.5518644 

0.5021333 

n.4632864 

0. 4323696 

0.407 3669 

0.3868654 

0.3698514 

0.3555831 

0.343 5072 

0.3332043 

0.3243518 

0.3166982 

0.3100448 


1.0000000 
0.8181402 
0.6924694 
0.6020795 
0.5349040 
0.4836044 
. n 4 :3 r -9! 
0.411 5170 
0.3856908 
0.3644028 
0.3466707 
0.3317335 
0.3190253 
0.3081179 
0.298 6S2 5 
0.2904632 
0.2832582 


1.0000000 
0.928253 3 
0.882489 1 
0.8529418 
0.8339617 
0.8221112 
VI '209° 
0.81i 8257 
0.8109901 
0.8120310 
0.8144719 
0.8179682 
0.8222664 
0.8271774 
0.8325582 
0.8382994 
0.8443164 
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